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MOTIVATIONS: The optimization of a structure and its control system has 

traditionally proceeded along two separate but sequential paths. First, the 
structure is optimized by selecting a set of member sizes a* which minimize a 
structural criterion J (a), subject to constraints h„(a), Eq (1). Then having 
specified the optimal structure, one may use the control theory to determine an 
optimal set of control variables if* that optimize a control criterion J c subject 
to constraints h c (uK Eq. (2). This two-step optimization procedure is the 
so-called separate optimization and is equivalent to finding the linear sum of two 
separate minima, Eq. (3). The question arises then as to whether it is possible to 
achieve a superior combined optimum (a**, u**) over (a*, u*) had one combined the 
two problems before, Eq. (4), rather than after, Eq. (3), the minimization. 
Intuitively, the answer to this question is affirmative, Eq. (5), since the 
minimum of the sum is less than the sum of the minima. 



• SIMULTANEOUS OPTIMIZATION 

J(a**, u**) = min [J (a) + J c (a, u)] (4) 

a, u 5 u 

J(a“, u**) < J (§*, y*) (5) 
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DEFINITIONS: A common starting point for most approaches begins with the 

second order dynamical equation, Eq. (6), in n $ degrees of freedom (d.o.f.). The 

M, D, K matrices are the mass, damping, and stiffness. G 1 = the disturbance 

influence matrix, G 0 = control influence matrix, a » structure design variables, u 
^-2 

= control variables, w = disturbance vector, and v = physical d.o.f. 

Let x = (v, v) T - state variables, then the equation of state is given by 
(7), with output consisting of controlled states z and measured states y, both 
related to x by Eq. (8). Where A, B, , B 0 , C, and C 0 are defined by Eq. (9). 


M(a) v+ D(a) v* + K(a) v = G ^w + G2 u 

i 

x = (v,v) T 

x = A(a) x + B 1 (a) w + §2 (a) u 

CONTROLLED STATES z = C 1 x 
MEASURED STATES y = C 2 x 


A = 


0 1 

1 

1 

1 

L 

- M ’ 1 D 



0 

-M ' 1 G 1 



Ci = (Cn,Ci2) ; O 2 =(C21.C22) 


( 6 ) 


(7) 

( 8 ) 


(9) 
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FORMULATION: Herein, we focus on LQ-based formulation as a natural one to 

generalize to the simultaneous control-structure optimization. Two types of 
controllers are considered; state feedback and output feedback. For both of these, 
the control criterion is taken as a quadratic function of the structural 
response and control energy, Eq. (10). For the structural criterion, we assume one 

that depends only on the structural vari abl es au As will become clear later, thi s 

simplifies the derivations considerably. An example of such structural criteria is 
the mass of the structure, M(a). 

The simultaneous optimization problem consists of finding the structure and 
control variables (a**, u**) that minimize the combined criterion (11), subject to 
any behavioral constraints (12), and/or side constraints (13) providing upper and/or 
lower bounds on the design variables (a, ui). Since the terms in (11) do not have 
the same units, the scalar a and matrices Q and R can be chosen on computational 
and physical grounds. 


• STATE FEEDBACK 

• OUTPUT FEEDBACK 

• INDIVIDUAL CRITERIA: 

CONTROL: J(a,u)=j (x T Qx + u T Ru)dt 

c~~ 0 ~ 

( 10 ) 

STRUCTURE: J s (a)=M(a) 

• COMBINED CRITERION: FIND (§**, y.**), a **e§, u“ e u 

J(a,u) = min[aM + | (x T Qx + u T Ru)dt] (11) 

a, u 0 ~ ~ ~ ~ ~ 

SUBJECT TO: hj (a, u)>o (12) 

a<a.<a, u<u.<G (13) 
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STATE FEEDBACK: With the assumption that the structural objective M is 

dependent upon a only, Eq. (11) simplifies to (14) and then to (15). This allows 
the familiar analytical solution for the optimal control u** in (18) and its 
companion equations (16) and (17). The necessary conditions for the minimum of (15) 
subject to the constraint (17) and the constraints imposed by (12) can be derived by 
first forming an auxiliary Lagrangian function, then setting its partial derivatives 
at the local minimum to zero. This yields conditions (19), (20) and (21), from 
which the optimal a** can be computed iteratively. For a given a, Eqs. (21) and 
(22A) are solved for P and P, a . With these, Eq. (19) is solved for an updated a, 
and (20) is evaluated to check the constraints. 


• J(a,u) = min|^o(M(a) + nmn J(x T Qx + u T Ru)dt] 

(14) 

= min|"aM(a) +Tr(PQ 0 )] 

(15) 

q 0 =b 1 b{ IF distw IS UNIT IMPULSE 
=x 0 xj IF DIST IS INITIAL COND x(o) = x Q 

(16) 

P SATISFIES A T P + PA +Q-P(B 2 R _1 B 2 )P=o 

(17) 

u‘* = -R -1 B 2 PX 

(18) 


• CONDITIONS OF OPTIMALITY 


ia. = al ^’a. +Tr ( PQ o’a +P - a Q 0 ) + 1 1 i h j-a =° 
i i i i j J 1 i 

X.h. =0 
J J 

a t p-.pa+q-pb 2 r _, Bj P=o 

AjP.a + p 'a A C +[A, J P+PA, a -P(B FT’bL,., P+Q a ]-0 
i i i i i i 

a c =a-b 2 r _1 bJp 


(19) 

( 20 ) 

( 21 ) 

(22)A 
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OUTPUT FEEDBACK: To avoid the state reconstruction necessary to implement the 

full state feedback control design, the static output feedback approach requires 
only the output of the measured states y. With a controller of the form of Eq. 
(22B), in which the gain F is assumed to stabilize the structure, the combined 
criterion in (11) reduces to (23). Equation (23) is similar to its counterpart, Eq. 
(15) for the state feedback, except now the output feedback gain £ may be considered 
as an optimization variable in addition to a. Furthermore, P(a,F) - P T (a,F)>0 
satisfies the Lyapunov Eq. (24). 

Here again, the necessary conditions for the minimum of Eq. (23), subject to 
the constraints of Eqs. (24) and (12), can be found by forming the auxiliary 
Lagrangian function and setting its partial derivatives at the local minimum to 
zero. This leads to Eqs. (25) to (29), which must be solved iteratively. For a 
given (a., F), Eqs. (27), (28) and (29), respectively, allow the solution of £, the 
Lagrangian matrix multiplier £ and the scaler Lagrangians for each behavioral 
constraint hj. With these, Eqs. (25) and (26) yield an improved (a, F), and so on. 


• ASSUMES u = (F STABILIZES STRUCTURE) 

• FIND MIN. OF J(a, u) = min [« M + Tr(PQb)] 

a, F 

WHERE P(a, F) SATISFIES AjB + PA c +Q c =0 

6 c = A + B 2 F C 2 
Q c = Q + cJf t RFC 2 


CONDITIONS OF OPTIMALITY 

a M. a j + Tr[PQ 0 a .+ 2LPA ca .+ LQ c ,aj3 + Si j h j, aj = 0 
Tr[2LPA c F+ LQ c ,FJ + f j h j,F = ° 

aJp*pa c+ q c =o 


tAl + A c L + Q o -0 


ijhj.O 


( 22 ) B 

(23) 

(24) 


(25) 

(26) 

(27) 

(28) 
(29) 
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EXAMPLE 1: STATE FEEDBACK (Ref. 1) 

The cantilever beam shown is modeled by three finite elements with 
cross-sectional areas a = (a^ a 2 , a 3 ) T , and has six d.o.f. An initial 
deformation vector at the six d.o.f. x(0) = x Q is specified, and a control force u 
is applied at the free tip. The areas a and control u are to be determined so as to 
minimize Eq. (11) while maintaining a fundamental open-loop frequency «>0.10 
rad/sec. Rather than a first order minimization, it was found necessary for faster 
convergence to use a second order scheme based on modified Newton-Raphson 
iterations. For this purpose, the design variables a and multiplier 2 are 
obtained iteratively from the recursive relations in Eq. (30). 



x = (0.01 1 , 0.001 35, 0.037, 0.002, 0.0688, 0.0021 6) T 
co >0.10 rad/sec. (i,e. h = co 2 -(0.10) 2 >o) 

E = 9.56 10 1 °N/m 2 , p =1660 Kg/m 3 
DAMPING = 0.5% CRITICAL 


a 

= , 

a 

- -S 

J,aa + Zijhi,aa 

" i 1 1 " 

h j,a 

i 

r+ 1 


r 

h . 

j.a 

0 


J, a + ^ihi 


J "j & 


h : 


(30) 
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ALGORITHM: The algorithm begins with a feasible initial design and a step 

length s with which a line search in the direction of negative gradient is 
performed. This is continued until the minimum is reached, or until the constraint 
is violated. If the latter occurs, constrained minimization is employed with an 
initial estimate of the multiplier ft from 

ft = -CH T H] _1 H T J,a 

where H = [h, . ... h„ 1 

With this Eq. (30) is used. The constraint is checked continually. If the 
design moves away from the constraint, unconstrained minimization is reverted to. 
Thus the minimization process alternates between iterations which involve 
unconstrained minimization and iterations which involve constrained minimization as 
outlined below. 


EXAMPLE 2: OUTPUT FEEDBACK (Ref. 2) 

In this example, an active disturbance force is applied to the free end of the 

beam in the figure below. At the other end the beam is pinned and a control torque 

is applied there. The measurements consist of angular deformation and angular 

velocity at the free end. The design variables for minimization are the 

T T 

cross-sectional areas a = (a 1 , a 2 , a 3 ) and the gains F - (F ] , . . . F 4 ) 1 . No 
behavioral constraints are imposed. Other parameters of the problem are listed 
below. Since there are no constraints, the minimization algorithm is essentially 
similar to the unconstrained gradient search portion of the algorithm described 
previously. 


◄ 1m » 

◄ 1m ► 

F4 1m N 

? 

a i 

a 2 

a 3 


y OJ y 2 y 3 Q y 4 


EXAMPLE 2 - RESULTS: The numerical results in the table below compare a. and ^F, 

and the resultant mass and combined index J for the initial design and optimized 
design. A factor of three reduction in J is realized as a consequence of 
simultaneous optimization over a. and £. In the accompanying plots, the transfer 
functions of the initial and optimum design from disturbance to the controlled 
output show three orders of magnitude reduction in response. 


SIMULTANEOUS OPTIMIZATION 



INITIAL DESIGN 

SIMULTANEOUS OPTIMIZATION 

AREAS 

a i 

a 2 

a 3 


a i 

a 2 

a 3 


0.1 

0.1 

0.1 


0.02308 

0.01654 

0.03572 

GAINS 

F 1 

F 2 

F 3 

F 4 

F i 

F 2 

F 3 F 4 


-1.0 

-1.0 

0.0 

0.0 

-0.0306 

-1.2872 

-1,2387 +0.1 136 

MASS 

30.0 




7.4 



A 

3 min 

3.06 x 10 4 



9.20 x 10 3 





FREQUENCY RESPONSE 



FREQUENCY RESPONSE OF (a) INITIAL DESIGN AND (B) OPTIMUM DESIGN 
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CONCLUDING REMARKS: To show the feasibility of simultaneous optimization as 

design procedure, we have used low-order problems in conjunction with simple 
control formulations. The numerical results indicate that simultaneous optimization 
is not only feasible - but also advantageous. Such advantages come at the expense 
of introducing complexities beyond those encountered in structure optimization 
alone, or control optimization alone. Examples include: larger design parameter 

space, optimization may combine continuous and combinatoric variables, and the 
combined objective function may be nonconvex. 

Future extensions to include large order problems, more complex objective 
functions and constraints, and more sophisticated control formulations will require 
further research to ensure that the additional complexities do not outweigh the 
advantages of simultaneous optimization. Some areas requiring more efficient tools 
than currently available include: multiobjective criteria and nonconvex 

optimization. We also need to develop efficient techniques to deal with 
optimization over combinatoric and continuous variables, and with truncation issues 
for structure and control parameters of both the model space as well as the design 
space. 


• SIMPLE FORMULATIONS USED WITH LOW-ORDER PROBLEMS 


• RESULTS SHOW SIMULTANEOUS OPTIMIZATION FEASIBLE AND ADVANTAGEOUS 

• ADDITIONAL COMPLEXITIES: - LARGER PARAMETER SPACE 


- POSSIBLE NONCONVEXITY OF OBJECTIVE FN. 

- MIXTURE OF CONTINUOUS AND COMBINATORIC 
VARIABLES 


• FURTHER EXTENSIONS: 


LARGER PROBLEMS 


- OTHER OBJECTIVE FNs, CONSTRAINTS, MORE 
SOPHISTICATED CONTROL FORMULATIONS 


- MORE EFFICIENT TOOLS TO DEAL WITH ABOVE 
COMPLEXITIES 
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UNIFIED TRUNCATION METHODOLOGY FOR 
CONTROL & STRUCTURE PARAMETERS OF 
MODEL SPACE & DESIGN SPACE 
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